Electron correlations in metal nanoparticles with spin-orbit scattering 
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The combined effect of electron-electron interactions and spin-orbit scattering in metal nanoparti- 
cles can be studied by measuring splitting of electron levels in magnetic field {g factors) in tunneling 
spectroscopy experiments. Using random matrix theory to describe the single-electron states in the 
metal particle, we find that even a relatively small electron-electron interaction strength (ratio of 
exchange constant J and mean level spacing 5 ~ 0.3) significantly increases g-factor fluctuations for 
not-too-strong spin-orbit scattering rates (spin-orbit time Tso ^ 1/^)- In particular, factors larger 
than 2 could be observed. (This is a manifestation of the many-body correlation effects in nanopar- 
ticles). While so far measurements only on noble metal (Cu, Ag, Au) and Al samples have been 
done for which the effects of electron-electron interactions are negligible, we discuss the possibility 
of observing interaction effects in nanoparticles made of other metals. 
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I. INTRODUCTION 

While the study of the combined effect of electron- 
electron interactions and elastic impurity scattering in 
two dimensions and near the metal-insulator transition 
in three dimensions remains one of the most impor- 
tant problems in solid state physics, the description of 
electron-electron interactions in disordered normal metal 
nanoparticles {i.e., "zero dimensions") has been found to 
be remarkably simple^'^: 

At a fixed number of electrons and without spin-orbit 
scattering, the only relevant interaction term is a long- 
range exchange interaction'^ 



Hi, 



(1) 



that couples to the total spin S of the nanoparticle. The 
exchange constant J is closely related to one of the Fermi 
Liquid constants of the bulk metal, and is independent 
of the details of the impurity configuration inside the 
nanoparticle. The interaction Hamiltonian (1), which is 
known as "universal interaction Hamiltonian" , is the only 
form of the electron-electron interaction compatible with 
random matrix theory.^ Random matrix theory provides 
a valid description of single-electron states as long as the 
dimensionless conductance g of the nanoparticle, which 
is the ratio of the Thouless energy i^xh and the mean 
level spacing S, is large. "'■^ Residual interaction terms not 
included in Eq. (1) are sample specific and small in com- 
parison to Eq. (1) by, at least, a factor 1/g. 

In the presence of spin-orbit scattering, spin is also 
randomized, giving rise to both sample-to-sample fluctu- 
ations of the electron-electron interaction and a suppres- 
sion of the exchange interaction (1). Since the spin-orbit 
scattering rate 7so = h/Tso plays the role of a "Thou- 
less energy" for the spin degree of freedom, for strong 
spin-orbit scattering the residual exchange interaction 
becomes small by a factor ^so/S 3> 1 in comparison to the 
interaction strength without spin-orbit scattering if spin- 
orbit scattering. However, the exchange interaction re- 
mains the dominant contribution to the electron-electron 
interaction as long as 7so/<5 <C g. 

In this paper, we present a detailed analysis of the com- 



bined effect of spin-orbit scattering and electron-electron 
interactions in the regime of moderate spin-orbit scat- 
tering, ~ S. In this parameter regime, the exchange 
interaction is not fully suppressed, while fluctuations are 
of the same order as the average.® This makes the regime 
of moderate spin-orbit scattering rates qualitatively dif- 
ferent from that without spin-orbit scattering (730 = 0) 
and that of strong spin-orbit scattering {jso/S 3> 1.) 
The parameter regime 7so ~ (5 is of interest for recent 
experiments on metal nanoparticles, in which the 
magnetic-field dependences of many-electron levels has 
been measured using tunneling spectroscopy.^^ More- 
over, an analysis of the regime 7so/<5 ~ 1 serves as a 
model study of the breakdown of the "universal interac- 
tion Hamiltonian" when the dimensionless conductance 
is small. A brief account of some of our findings was 
previously published in Ref. 12. 

Experimentally, the spin structure of electronic states 
can be measured through the magnetic-field dependence 
of steps in the current-voltage characteristic of a metal 
nanoparticle coupled to source and drain electrodes via 
tunneling contacts. These steps occur if the applied volt- 
age is equal to the difference of the energies of many- 
electron levels \Ne + 1; k) and | A^e; I) which have A''e and 
A^e -I- 1 electrons, respectively,^^ 



(2) 



The derivative dVki/dB of the voltage at which a step 
occurs to the magnetic field B is parameterized through 
a "g-factor". 



dVki 1 



(3) 



where fig — \e\h/2mc is the Bohr magneton. Without 
spin-orbit scattering, but with electron-electron interac- 
tions, many-electron states are characterized by their to- 
tal spin S and by its z component. Since tunneling spec- 
troscopy measures transitions in which the electron num- 
ber changes by one, the total spin S of the nanoparticle 
changes by 1/2 upon addition or removal of an electron. 
This "selection rule" renders all observed g factors equal 
to two, irrespective of the spin Sk and Si of the two 
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many-electron levels participating in the transition. On 
the other hand, without interactions, the tunneling spec- 
troscopy g factors correspond to the magnetic moment 
of one single-electron level only, giving rise to a distribu- 
tion of g factors in which all levels have, at most, spin 
1/2.^'''^^'-'^^ As was shown in Ref. 12, the combined ef- 
fect of spin-orbit scattering and electron-electron interac- 
tions is to simultaneously suppress the spin of the many- 
electron states and lift the "selection rules", causing a 
much wider distribution of tunneling spectroscopy g fac- 
tors than in the non-interacting case. In particular, there 
is a significant probability to find g factors larger than 
two if spin-orbit scattering is not too strong (7so/^ ^ 2). 
The occurrence of g factors larger than two is a unique 
signature of the interplay of electron-electron interactions 
and spin-orbit scattering. 

What are the main differences between g factor dis- 
tributions with and without electron-electron interac- 
tions? In order to answer that question, we note that 
the electron-electron interactions has two main effects in 
a nanoparticle: to organize the many-electron states ac- 
cording to their total spin S, i.e., to lift the degener- 
acy between states of different S but with the same or- 
bital content, and to lower the energy of a many electron 
state with spin 5 by the amount JS{S + 1), increasing 
the abundance of high-S' states among low-energy ex- 
cited states, see Fig. 1. (In fact, for J/6 > 0.3, there 
is a significant probability that the ground state has a 
nontrivial spin S > 1/2, see, e.g., Refs. 17,18). It is be- 
cause of the combination of these two effects, together 
with the; lifting of selection rules by spin-orbit scattering, 
that electron-electron interactions enhance the width of 
the 5(-factor distribution so significantly. Moreover, be- 
cause the relative abundance of high-S* states depends 
on the excitation energy, the g factor distribution will 
be different for transitions to an excited state than for 
transitions to the groimd states. Again, this is different 
from the non-interacting case, where g factor distribu- 
tions for transitions to the ground state and to excited 
states are eqiial. A third difference between the cases 
with and without interactions is that, as the selection 
rules are gradually broken down by spin-orbit scattering, 
different transitions may have very different weights, in 
contrast to the non-interacting case, for which all tran- 
sitions have weights within a factor of order unity from 
each other. (The "weight" of the transition is the height 
of the corresponding step in the current- voltage charac- 
teristics.) 

In principle, one should consider contributions to the 
g factor from the orbital magnetic moments of the en- 
ergy levels and from the spin magnetic moment. In this 
work, we consider the spin contribution to the g factor 
only and neglect the orbital contribution. For the spin- 
orbit scattering rates jso/6 ^ 1 we consider here, this is 
justified if the electron motion in the metal nanoparticle 
is diffusive with mean free path I much smaller than the 
nanoparticle size L.^^'^^ 

The outline of this paper is as follows. In section 
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FIG. 1: Occupation of single-electron levels for the lowest- 
lying many-electron states with total spin 5 = (a) and 5=1 
(b), and for S = 1/2 (c) and 5 = 3/2 (d). In the absence of 
exchange interactions, the states with spin S — and 5 = 1/2 
arc the ground states for even and odd numbers of electrons, 
respectively. The exchange interaction compensates (part of) 
the kinetic energy cost of the higher spin states. 



II we present the theoretical formalism. Since we only 
consider spin-orbit scattering rates jso/S <C g, random- 
matrix theory can be used to describe the single-electron 
states. In Sec. Ill we discuss the results of numerical 
simulations for the g factor distribution for transitions 
from the iVe-electron ground state to the (iVe-l-l)-electron 
ground state and (A^e + l)-electron excited states, where 
Ne is taken even. The restriction to transitions start- 
ing from the A^g-electron ground state is appropriate if 
the metal nanoparticle relaxes to the A^e-electron ground 
state between tunneling events. In Sec. IV we discuss 
the consequences of our findings for various metals: our 
results depend on the ratio J/S, which strongly depends 
on the metal under consideration. We conclude in Sec. 
V. Finally, in the appendix, we report an analytical cal- 
culation of the g factor distribution for weak spin-orbit 
scattering, 7so/^ -C 1, again paying special attention to 
differences between transitions involving ground states 
only and transitions to or from an excited state. 



II. THEORETICAL DESCRIPTION 

A. Effective Hamiltonian 

Random matrix theory can be used to describe the 
single-electron wavefunctions and energy levels of a metal 
nanoparticle. For a metal nanoparticle with spin orbit 
scattering, the appropriate random matrix ensemble in- 
terpolates between the Gaussian Orthogonal Ensemble 
(GOE) and the Gaussian Symplectic Ensemble (GSE) of 
random matrix theory. 



Ho{X) = Hqoe + Hso{X). 



(4) 
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Writing the spin degrees of freedom explicitly, one has 

-f^GOE = 5012, (5) 

iX ^ 

where I2 is the 2x2 unit matrix in spin space, <jj is 
the Pauli matrix (j = 1,2,3), S is a,n N x N real sym- 
metric matrices, and Aj is a real antisymmetric matrix 
{j = 1,2,3). The elements of the matrices S, Ai, A2, 
and As are drawn from independent Gaussian distribu- 
tions with zero mean and with equal variances for the 
off-diagonal elements. The diagonal elements of S have 
double variance, whereas the diagonal elements of Ai, 
A2, and A3 arc zero because of the antisymmetry con- 
straint. In random matrix theory, the limit ^ 00 is 
taken at the end of the calculation. The random ma- 
trix description is valid as long as energy differences of 
the energy levels and wavefunctions of interest are small 
compared to the Thouless energy i^xh-'*'^ For a disor- 
dered metal nanoparticle of size L, mean free path /, and 
Fermi velocity vp, Erh ~ vpl/L"^- 

The parameter A in Eq. (4) determines the strength of 
the spin-orbit scattering, = n'yso/S = tt/tsqS, where 
Tso = l/7so is the spin-orbit time and S is the mean spac- 
ing between (spin-degenerate) eigenvalues of -ffcoE, sec 
Refs.^*'"'^^'^^ The case A = corresponds to the absence 
of spin-orbit scattering, while the limit A — > 00 describes 
the situation when spin-rotation symmetry is completely 
broken. The factor l/VlV in front of Hso ensures that 
the relation between A and the physical parameters Tso 
and S does not involve the matrix size A''. 

Each eigenvalue s:^ of the Hamiltonian (4) is doubly 
degenerate (Kramers degeneracy). After diagonalization, 
the Hamiltonian Hq can be written as 

/?0=E^/'(^ilV+V'i2V'^2), (7) 

where ipjj^^ and -0^^ are creation and annihilation oper- 
ators for an electron in the state \fia), where a = 1,2 
labels the two time-reversed states in the Kramers dou- 
blet. 

Combining the single-electron Hamiltonian (7) and the 
interaction Hamiltonian (1), one find the total Hamilto- 
nian 

^ = E ^MlA^ + - js^- (8) 

Equation (8) is valid up to a charging energy that de- 
pends on the electron number Ne only; the charging en- 
ergy plays no role for in the problem we consider. 

In the absence of spin-orbit interaction, the exchange 
interaction (second term in Eq. (8)) commutes with the 
non-interacting Hamiltonian Hq (first term in Eq. (8)). 
The many-electron eigenstates are found by diagonalizing 
Hq at a fixed value of the total spin S and its z compo- 
nent Sz- With spin-orbit interaction, however, the inter- 
action does not commute with Hq. This has important 



consequences for the ground state and for the excitation 
spectrum of a metal nanoparticle. Typically, for most 
normal-mctal nanoparticles and for quantum dots, J is 
estimated to be in the range < J/(5 < 1, see Sec. IV. 
While this implies that the effect of the exchange inter- 
action cannot be treated in first-order perturbation the- 
ory, interaction effects only cause correlations in a small 
window around the Fermi energy which is, in principle, 
available to direct numerical diagonalization. Hereto, we 
write the operator S for the total electron spin in terms 
of the creation and annihilation operators ijjj^^ and i^^^a 
of the single-electron Hamiltonian Hq, 

S = E E V'iaS/^a,L'/3V'!^/3, (9) 
IJ,,v a, 13=1,2 

where 

{si)na,uf3 = '^Wf3\a-i\iia), i = 1,2,3. (10) 

The quantity of interest in our calculation is the 
magnetic-field dependence of the many-electron energy 
levels for small magnetic fields and an odd number of 
electrons, which is described through the g factors, see 
Eqs. (2) and (3) above. The magnetic-field dependence 
arises both through the Zecman coupling to the elec- 
tron spin and through the orbital coupling to the angu- 
lar momentum. ^^'^^ For large diffusive metal nanoparti- 
cles and for not-too-large spin-orbit scattering strengths 
TsoS > 1, the Zeeman coupling dominates. ^^'^^ Since the 
interaction effects studied here are most important for 
TsoS ^ 1 (sec below), we neglect the orbital contribution 
to the g factors in the discussion below. For a magnetic 
field H along the z axis, the Zeeman coupling to the 
magnetic field is described by the Hamiltonian 

Hz = -2iiBHSs, (11) 

where the ^;-component of the total spin is given by 
Eq. (9) above. 

B. Tunneling spectroscopy 

Following Ref. 13, we assume that the conductance of 

the tunneling contact connecting the nanoparticle to the 
source reservoirs is much smaller than the conductance 
of the contacts connecting the particle to the drain reser- 
voir, so that the current / through the nanoparticle is lim- 
ited by the processes where an electron tunnels onto the 
particle. In this case, one can assume that relaxation is 
sufficiently fast that the nanoparticle is in the Afg-particle 
ground state \Ne, 0) before each tunneling event. It fol- 
lows that current steps occur in the current- voltage char- 
acteristic or, equivalently, a peak in the nanoparticle's 
differential conductance dl/dV, when the source-drain 
voltage eV = eVko = i?Ar^+i,fe — En^^, sec Eq. (2). For a 
point contact that injects electrons into the nanoparticle 
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H=0 



(e) 



FIG. 2: Top panel: schematic representation of four lowest 
many-electron levels for a nanoparticle with an odd number of 
particles without spin-orbit scattering. The spacings between 
the single-electron levels are as indicated in the figure. The 
figure represents the case < Aq -|- Ai — 3J < Ao, Ai. Note 
that the first excited state has spin 5* = 3/2 and is fourfold 
degenerate. Bottom panel: Evolution of energy levels of the 
nanoparticle for small spin-orbit scattering rates A and mag- 
netic fields H. Spin-orbit scattering separates the quadruplet 
(b) into two doublets with ill-defined spin. The magnetic field 
lifts the degeneracy of all remaining doublets. 



at position r, the size of the current step is proportional 
to the matrix element 



\{Ne 



l,k\^]{r)\N„0)f 



|(7V, + l,fc|^|(r)|iV„0)p, 



where the creation operator ipl-{r) creates an electron 
with spin a in the nanoparticle at position r. In terms of 
the basis of single-electron states, one has 



(13) 



In random matrix theory, the matrix element {ra\fia), 
a = ±1, is replaced by an (arbitrary) spinor in the 2N- 
component vector representing the state fia, a = 1,2. 

In the presence of both spin-orbit scattering and 
electron-electron interactions, the iVe-electron states 
\Nf,,l) are non-degenerate for zero magnetic field if Ng 
is even. In that case, dEi/dH = at iJ = 0. On the 
other hand, (iVg -I- l)-electron states are doubly degener- 
ate if Nf, is even. These states split in a magnetic field. 
For small magnetic fields the splitting is linear. Hence, 
if TVe is even, the g factor g^i associated with voltage V^i 
at which a current steps occurs is directly related to the 



magnetic-field derivative of the {N^ 

|7Ve + l,0, 



dE 



dB 



1 



No, 



l)-electron level 



(14) 



In the remainder of this paper we continue to refer to the 
even-electron ground state as the "iVe-electron ground 
state" and to the odd-electron states as the "(iVg + 1)- 
electron states" . 

An example showing how the various degeneracies are 
lifted by the exchange interaction and by a magnetic field 
is shown in Fig. 2. The top panel of the figure shows the 
four lowest many-electron states for an odd number of 
electrons, for the specific case that the ground state has 
spin 5 = 1/2 and the lowest excited state has 5 = 3/2 
in the absence of spin-orbit scattering. Without spin- 
orbit scattering, the 5 = 3/2 state is fourfold degenerate. 
Spin-orbit scattering lifts the fourfold degeneracy of the 
5 = 3/2 quadruplet, separating it into two doublets with 
ill-defined spin (lower panel of Fig. 2, center). Finally, an 
applied magnetic field lifts the degeneracy of all doublets. 
For an even electron number, spin-orbit scattering lifts 
all possible degeneracies; to first order in the field, the 
applied magnetic field has no effect. 

The definition that g factors are derivatives of energy 
levels to the magnetic field implies that the Zeeman en- 
ergy scale is the smallest nonzero energy scale in the 
problem. In particular, it is smaller than the spin-orbit 
induced splittings of high-spin states (see Fig. 2), these 
splittings being proportional to the spin-orbit scattering 
rate oc I/tsq. However, when g factors are calculated 
without spin-orbit scattering, the Zeeman energy is (by 
definition) larger than the spin-orbit rate. We'll find be- 
low that the two limits do not commute in the presence 
of electron-electron interactions, i.e., that g factor distri- 
butions calculated in the limit A — *■ 0-|- are different from 
the g factors at A = 0. (Without spin-orbit scattering all 
peaks in the differential conductance split with g factor 
g = 2.) It should be pointed out that experiment involves 
finite magnetic fields, for which the Zeeman energy can 
be larger than the spin orbit rate. For such magnetic- 
field dependences of the levels, a different slope at zero 
field can easily go unnoticed. 

Although g factors contain information on the posi- 
tions of peaks in the tunneling spectrum of the nanopar- 
ticle, knowledge of the sizes of the peaks in the differential 
conductance is important for a correct interpretation of 
the results. Peak heights contain information that would 
be formulated in terms of "selection rules" in the ab- 
sence of spin-orbit scattering. Indeed, without spin-orbit 
scattering, the total spin 5 and its z component 53 can 
change by 1/2 at most in a tunneling process, which lim- 
its the possible transitions between and + 1 elec- 
tron states, and, hence, the possible locations of peaks in 
the differential conductance. It is because of the selection 
rules that one does not observe g factors larger than two, 
despite the fact that there exist high-spin many-electron 
states. Similarly, without interactions, the occupation 
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selection rules. The figure shows two possible A'e-electron 
ground states (left) and three [N^ + l)-particle excited 
states (right). Without spin-orbit scattering or without 
exchange interactions, the states (c) and (d) can be ac- 
cessed from the A^e-particle ground state (a), but not 
state (e). Similarly, state (e) can be accessed from ground 
state (b), but not (c) and (d). When spin-orbit scatter- 
ing and exchange interactions are both present, the A^e- 
elcctron ground state is a superposition of the states (a) 
and (b) and all possible transitions have a finite matrix 
element. Notice that, since the energy difference between 
states (a) and (b) is typically small — on average, equal 
to (5 — 2 J, mixing of these two states is strong for spin- 
orbit scattering rates A 1. 



FIG. 3: Loft: Sclicmatic representation of ground states for 
an even number of electrons without spin-orbit scattering. 
Depending on tlie strength of the exchange interaction, the 
A^e-electron ground state may have spin 5 = (a) or S = 1 
(b). Right: three excited (ATg -|- l)-electron states. 



of single-electron levels cannot change by more than one 
electron, which also limits the number of allowed transi- 
tions in tunneling spectroscopy. With spin-orbit scatter- 
ing and interactions, peaks that were previously "forbid- 
den" are present, in principle, although their height may 
be small. We return to this issue in more detail in Sec. 
III. 

The example in Fig. 3 may further clarify the role of 



C. Matrix elements of interaction Hamiltonian 

For the construction of the interaction Hamiltonian in 
the basis of many-electron eigcnstates of the Hamiltonian 
Hq, one needs explicit equations for the matrix elements 
of the exchange interaction in that basis. Since the ex- 
change interaction changes the single-electron states of at 
most two electrons, the only nonzero matrix elements of 
the exchange interaction occur between states that can be 
written in the form i>l,ai'ij3\F) and V'jj'a'V'I'/S'l-^)' where 
\F) a certain reference non-interacting state. After some 
algebra one then finds 



{F\i^.>04,'c.'S^lJl0\F) 

(15) 



where 



(F|S^|F) = -(A^e — 2) — 2 ^ ^ [s07,,^'7' • s,p'j>^^^ — s^-y,^^ • s^'Y^^'ji] . (16) 

0,7eF 4>',j'eF 



I 

The summations over (j) and 7 extend over all single- ^jjial^) and ^/j^^jF). For those states, one needs the ma- 
electron states 1^,7), 7 = 1,2, that are occupied in the trix elements 
state \F). 

Similarly, matrix elements of the Zeeman energy Hz 
are nonzero only if many-electron states differ by not {F\7p^i3Si)lJF) = s^p^^a + 5^0,^^ ^07.07- (1^) 

more than one electron, i.e., between states of the form <P,jeF 
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At this point, it is important to verify the apphcabihty 
of the random matrix theory. In order for random matrix 
to apply, summations over the Fermi sea should converge 
within a Thouless energy from the Fermi level. In Eqs. 
(15), (16) and (17), the sums of the form J2,j> jeF^'Pi^'Pi 
clearly satisfy this condition, by virtue of the equality 

S01,01 + S02,02 = 0, (18) 

which follows from the observation that the states |(/), 1) 
and 10,2) are time reversed. The sums of the form 
^0 ^^p Si^a_4>-f ■ ^4>-i,vi3 also meet this condition, since the 
summand decreases oc l/{eF — £4,)"^ for </> well below the 
Fermi level and n and v close to the Fermi level. In 
the diagonal matrix element (16) the sum over (</',7) 
and (0',7') is logarithmically divergent as a function of 
the Fermi energy. This, however, has no consequences 
for the magnetic-field dependence of the many-electron 
states and the peak heights, since the divergence is the 
for all matrix elements and simply corresponds to the 
overall shift of the ground state energy. We conclude that 
random matrix theory can be used to access the many- 
electron ground state and the low- lying excited states. 



between tunneling events, the (statistical) properties of 
the A^g-particle ground state play a key role in determin- 
ing the g factor distribution. Without spin-orbit scatter- 
ing for J < 0.3(5, there is only a small probability that 
the A^e-particle ground state has spin one, and a vanish- 
ing probability that the A^e + 1-particle ground state has 
spin 3/2 (Ne is assumed even).^'''^^ The probability to 
find an A^e-particlc ground state with 5 = 1 becomes ap- 
preciable for J > 0.3, whereas the probability to find an 
(iVe + l)-particle ground with spin 3/2 becomes signifi- 
cant for J > 0.55 only. For the values of J we consider, 
states with spin > 5/2 do not play a role; they have been 
excluded from the truncated many-electron basis. 

The strength of spin-orbit scattering strength A is 
taken from to 2.8. Although larger spin-orbit scat- 
tering strengths do occur in metal nanoparticles,^'^" in- 
teraction effects are small at those values of A and the 
non-interacting theory of Refs. 14,15,16 works well. 

The random matrices in our simulation are taken of 
size 2A'' = 400. This ensures that the condition \^ <^ N 
necessary for the applicability of the random matrix (4) 
is satisfied for all values of A. For A < 2 we have taken 
2^" = 200 in the simulations. 



III. RESULTS AND DISCUSSION 

In order to calculate g factors and peak heights, we 
have diagonalized the Hamiltonian (8) numerically. 

Our numerical procedure is as follows: We first diago- 
nalized the non- interacting Hamiltonian Hq. The inter- 
action is considered in a truncated basis of the many- 
electron states, taking the 92 lowest lying many-electron 
eigenstates of Hq for a A^e + 1 electrons (A'e even), and 
the 76 lowest lying states for Ne electrons. Finally, we 
diagonalized the truncated interaction and found the g 
factors of the M = 8 lowest-lying (A''e + l)-electron states, 
and the peak heights that follow for transitions from the 
A'e-electron ground state. For the calculation of the g 
factors we introduce a small magnetic field and calculate 
the magnetic field derivative numerically. We verified 
that truncating the interaction Hamiltonian at the low- 
est lying 92 and 76 many-electron states has no effect on 
the final results by comparing our results to those that 
were obtained using a smaller basis set. 

We have investigated exchange-interaction strength 
ranging from J = to J = 0.6(5 which is valid for 
most metals, see Sec. IV. The same parameter range 
should apply to quantum dots. Analysis of the Coulomb 
blockade peak spacing distribution^^ suggests that J/5 
is between 0.3 and 0.4 in large quantum dots in a 
GaAs/GaAlAs heterostructure; whereas recent density 
functional studies of ground state spin distributions in 
ballistic quantum dots are compatible with the Hamilto- 
nian (8) only if J « 0.6(5, see Ref. 20. 

Important changes occur within the range of exchange 
interactions we address here. Since the metal nanopar- 
ticle is assumed to relax to the A^e-particle ground state 



A. Average g factors 

We have calculated ensemble averages of the g factors 
{gk), fc = 0, 1, . . . ,M-1 of the M lowest (Afe-hl)-electron 
states. Here gk is the g factor corresponding to the fcth 
(A^e + l)-electron system, fc = 0, 1,...,M-1. The en- 
semble average is taken over 300 realizations. In Fig. 4 
we show the ensemble averaged g factors for the ground 
state and the first excited state, {go) and (gi), as well as 
the average over all calculated g factors {g) 

M-l 

g = M-^Y. 9k- (19) 

fc=0 

The first observation to be made from Fig. 4 is that 
for J > 0.4(5 and A < 2 interactions lead to a significant 
increase in the average g factor. In fact, there is a sig- 
nificant parameter range for which {g) > 2. The origin 
of the large g factors is that exchange interactions lift 
the degeneracy with respect to the total spin S. Hence, 
with exchange interactions, many-electron states with a 
finite spin are energetically separated from states with 
spin or 1/2. For the parameter range considered here, 
the relevant nontrivial spin states have 5 = 3/2 for an 
odd number of electrons. The role of spin-orbit scatter- 
ing is to lift the fourfold degeneracy of the 5 = 3/2 states 
and, for larger spin-orbit strengths, to suppress the spin 
content of the single-electron states that build the many- 
electron state. Let us first discuss the effect of lifting the 
degeneracy of the 5 = 3/2 state by spin orbit scattering. 

In general, spin orbit scattering splits the fourfold de- 
generate 5 = 3/2 state into two doublets. Neglecting 
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contributions from other many-electron states, each dou- 
blet consist of two states that can written in the form 



3/2 



|state) = a„|3/2,n), 

ri=-3/2 



(20) 



and the time-reversed of Eq. (20), which is obtained by 
sending a„ sign(n)a!l„. Here IS*, 5*3) is the {Ng + 1)- 
electron state with total spin S and z component of the 
spin equal to 6*3. Because the spin-orbit matrix elements 
are random, the amplitudes a„ are essentially random as 
well. (This statement is verified in the appendix.) The g 
factor of the state (20) is 




3/2 

E 

Ti=-3/2 



4|n|a„a_ 



(21) 



One easily verifies that this can be larger than two. With 
exchange interaction but without spin-orbit scattering, 
there is a finite probability that the A^e-particle ground 
state (Ne even) has spin S = 1. In that case, it has two 
singly occupied orbitals and, hence, in principle, a finite 
overlap with a state of the form (20) after addition of 
an electron. With spin-orbit scattering, the TVe-electron 
gromid state is guaranteed to be non-degenerate, so that 
its derivative to the magnetic field is zero. We conclude 
that, the g factor of the state (20) can be larger that 
two, that it can correspond to a transition between the 
A^e-electron ground state and an (A^e + l)-electron state, 
and that the corresponding conductance peak has a finite 
height. 

A finite amount of spin-orbit scattering stabilizes the 
above arguments by increasing the splittings between 
many-electron states that are degenerate in the absence 
of spin-orbit coupling. On the other hand, with moder- 
ate spin-orbit scattering, more many-electron states are 
added in the doublet (20). This has two consequences: 
(1) the the spin content of each of the underlying single- 
electron states is reduced, which, eventually, leads to 
a suppression of g factors, and (2) when more many- 
electron states are admixed, overlaps and, hence, peak 
heights are increased, so that the role of selection rules is 
further diminished. In order to illustrate the value of A 
needed to admix different many-electron states, we note 
that for the A^e-particle ground state without spin-orbit 
interaction, the energy separation between the S = 1 and 
S = ground states is 2 J — 5 on average. Hence, even 
for a relatively small spin-orbit scattering rate A ~ 0.5 
the ground state with spin-orbit scattering will have sig- 
nificant weight in both of these states. 

Also notice that, unlike in the non-interacting case, the 
average g factor depends on the excitation energy of the 
many-electron state for A < 2, see Fig. 4. The origin 
of this dependence is that, without spin-orbit scatter- 
ing, the probability of that an {Ng + l)-particle state has 
nontrivial spin {S > 3/2) increases with the excitation 
energy. As discussed above, although spin-orbit scat- 
tering lifts the fourfold degeneracy of these states and 
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FIG. 4: Ensemble averaged g-factors. Top: average g factor 
(go) of (A'e -I- l)-particle ground state for J = (the lowest 
solid curve), J = 0.35, J = 0.45, J = 0.55, and J = 0.65 Mid- 
dle: average g factor {gi) of first excited (A^e-l-l)-particle state 
for J/5 = 0.1, 0.2, 0.3, 0.4, 0.5 and 0.6. Bottom: ensemble av- 
eraged g factor averaged over the first M — 8 (iVe-l-l)-electron 
states for J/5 = 0.1, 0.2, 0.3, 0.4, 0.5 and 0.6. 



suppresses the spin, it is the underlying nontrivial spin 
character persisting to finite A that gives rise to the in- 
creased g factors. 

A remarkable feature of Fig. 4 is that the ensemble 
averaged g factor (g) docs not approach two in the limit 
A ^ 0. On the other hand, without spin-orbit scattering, 
all observed g factors should equal g = 2. In Sec. II we 
discussed why there can be a difference between g fac- 
tors in the limit A ^ and g factors calculated without 
spin orbit scattering, i.e., at A = 0. In fact, in Fig. 4, 
(g) is overestimated for A ^ 0, because the plain ensem- 
ble average does not take tunneling spectroscopy peak 
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FIG. 5: Top: average of all calculated g factors, where each g 
factor is weighed by its normalized peak height (22). Bottom: 
average of all calculated g factors for which the normalized 
peak height is larger than 0.1. In both panels, results are 
shown for J/d = 0, 0.1, 0.2, 0.3, 0.4, 0.5, and 0.6. For J = 
there is no difference with Fig. 4. Inset of lower panel: the 
probability for a level to be visible in the experiment, i.e. to 
have a weight larger than the threshold one, see (23). 

heights or "selection rules" into account: The average 
is taken over all (Ne + l)-particle states irrespective of 
the height of the corresponding peak in the differential 
conductance. In particular for small J, one expects that 
(7Ve + l)-electron states with g factors larger than two are 
likely to have small tunneling spectroscopy peak heights. 

We have taken two different approaches in order to 
account for "selection rules" . First, we have replaced the 
ensemble average by a "weighed" average, in which every 
g factor is weighed by the normalized peak height 

\k=o I Z^fc=o Wk 

In the second approach, we have removed all peaks with 
normalized weight Wk below a certain threshold value, 
where we arbitrarily set the threshold to W\, = 0.1 x 
max^i;^ Wk, 

/ 1 if Wfc > Wt, . . 

Wk^Wk,^ = ^Q if«,fe<«;t. (2^) 

In this method, the number of levels per realization de- 



pends on the realization, 

M^ = Y,Wk,t (24) 

k 

and the average ^-factor is determined through 

I ^ M-l \ 

(5) = ( Ti- X] ) • (25) 

\ * fc=0 / 

The threshold mimics the experimental reality that small 
peaks cannot be distinguished from the noise, and, hence, 
have their g factors left out in the statistical analysis. 
While the second method is closer to the way experi- 
ments are analyzed - all g factors of conductance peaks 
that are observed are taken equally into account in the 
average — the first method has the advantage that it 
does not contain the somewhat arbitrary threshold at 
Wk = vJt = 0.1. Both methods enforce the "selection 
rules" in the absence of spin-orbit scattering. They also 
give almost identical results for the average g factor, as 
is seen from Fig. 5 where we show the weighted ensemble 
average of the g factors of all levels considered (g) , as well 
as the ensemble average calculated using the "threshold" 
method. As shown by comparison of Figs. 5 and 4, in 
the limit A ^ 0, the average g factors are close to two 
for small J, whereas {g) is significantly higher than two 
for J > 0.3. In the inset of the lower panel of Figs. 5 
we show the probability of the level to be visible, i.e. to 
have a peak height above the threshold. Remarkably, the 
curves for J > 0.35 have a maximum for moderate values 
of spin-orbit scattering A ~ 1. This has a direct physical 
interpretation: both for small and large A, approximate 
selection rules are in place. For small A these selection 
rules represent the conservation of spin at A = 0, whereas 
for large A they follow from the suppression of the (ex- 
change) interaction, which causes the remaining physics 
to be single-particle like. 

B. fluctuations of g factors 

Cumulative probability distributions of g factors are 
shown in Fig. 6 for J = 0, A = 0.7, for J = 0.35, 
A = 0.7 and for J = 0.3^, A = 0.9 (we have taken into 
accoimt only peaks whose weights are nonzero according 
to the criterion (23)). Comparing the two distributions 
at A = 0.7, one notices that the exchange interaction has 
little effect on the tail of the g factor distribution for very 
small g factors. However, for larger g factors, the weight 
of the probability distribution is shifted towards larger g 
factors, including a long tail in the region g > 2. Fig. 6 
confirms the previous observation that the effect of ex- 
change interactions is to increase the average g factors. 
The spin-orbit scattering rate for the third probability 
distribution shown in Fig. 6, A = 0.9, has been chosen 
such that the average g factor {g)t « 1.58 coincides with 
that of the case J = 0, A = 0.7. Comparing the two 
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FIG. 6: Cumulative probability distributions function for g- 
factors for the cases A = 0.7, J = (dashed curve) A = 0.7, 
J = 0.3 (solid), and A = 0.9. ,/ = 0.3 (dotted). Even a 
relatively small strength of the exchange interaction is enough 
in order to broaden the distribution function significantly. 

probability distributions, we conclude that the interac- 
tions still lead to a significant increase of the g factor 

fluctuations, including a large probability to find g fac- 
tors larger than two, even if the average is well below 
two. 

In Fig. 7 we show the probability for a level to have 
a y-factor larger than two. In an experiment, typically 
g factors of 5-10 consecutive levels can be measured. ^'^^ 
From Fig. 7 we then conclude that there is a significant 
probability that one of these g factors is larger than two 
if J > 0.2(5. The bottom panel shows the probability 
that a level |A;) has a g factor larger than two and a 
weight Wk > Wt- As a result of the breakdown of selection 
rules, this probability increases with increasing spin-orbit 
scattering in the region A < 0.5: States which have large 
5-factors but small weights for small A become visible for 
larger values of A. The ratio of the probability shown in 
the bottom plot of Fig. 7 and that in the top plot of Fig. 7 
is the probability that a random level can be resolved in 
the experiment, see Fig. 5 

In addition to addressing the full probability distribu- 
tion of g factors, we should consider the possibility of 
correlations between g factors within one realization. In 
principle, such correlations can exist, because, within one 
realization, all g factors correspond to transitions from 
the same A^e-electron ground state (iVg even). Although 
the A^e-electron ground state does not affect the values 
of all possible g factors for the (TVe + l)-particle levels, it 
docs affect the peak heights, and hence determines which 
g factors possibly "drown" in the noise. 

In order to quantify g factor correlations, we have 
looked at the correlation function 

C(J,A) = /j2^(g,-{-g))(g,-{g))\. (26) 

For the calculation of the correlation function C(J, A) we 
removed all (^-factors with normalized weight Wk below 
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FIG. 7: Top: Probability for the g factor of a level to be 
larger than 2.0 for the values of the exchange constant J = 
0.1 -^ 0.65 as a function of the spin-orbit scattering strength 
A for a random visible (i.e., satisfying the threshold criterion, 
see Sec. Ill A) level. Bottom: Probability for a random level 
to be visible and, at the same time, have g > 2. for the values 
of the exchange constant J = 0.1 -t- 0.6(5 as a function of the 
spin-orbit scattering strength A. 



the threshold value Wt = 0.1 from the average, which 
means that the number levels Mt considered in the sum- 
mation becomes dependent on the actual realization. For 
the range of exchange interactions J and spin-orbit scat- 
tering rates A we considered, the correlation function 
C( J, A) was nonzero, but always smaller than 0.1. The 
maximum value C( J, A) ~ 0.1 was obtained for A ^ 0.5. 
Comparing the above difference with the typical variance 
of £f-factors, see Fig. 6, we conclude that correlations be- 
tween different levels within the same nanoparticle are 
not important if the number of g factors measured in a 
single metal nanoparticle does not exceed 10. 



C. Probability of nontrivial ground state 

In Fig. 8 we show the probability that the metal 
nanoparticle is found in the non-interacting ground state 
(Fermi sea), as a function of A and J. For A^e-electron 
states, the probability to find the nanoparticle in the non- 
interacting ground state deviates quite significantly from 
1 if J > 0.4 and A = 0.^^'^^ Upon increasing A, the prob- 
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see Ref. 21. The Fermi liquid parameter Fq appears in 
the expression for the paramagnetic susceptibihty 



FIG. 8: Probability to be in the non-interacting ground state 
for J = 0.3 — 0.65. The horizontal line P = 1 corresponds to 
the case J = 0. 



ability to be in the non-interacting ground state increases 
and approaches unity when the spin-orbit scattering rate 
exceeds the exchange interaction. 



IV. MATERIAL DEPENDENCE 

Petta and Ralph have measured the probability distri- 
bution of g-factors of Cu, Au, and Ag nanoparticles.^'^° 
Theoretical estimates and experimental investigations of 
the exchange interaction in the noble metals show that 
J/6 < 0.1 for Cu, Au, and Ag.^^'^^ Hence, the inter- 
action effects in the above materials are very small and 
it is natural that the existing experiments can be ex- 
plained quantitatively using theory for the noninteract- 
ing (J = 0) case. Indeed, both the average and the width 
of the g factor probability distribution measured in Rcfs. 
9,10 were found to be in good agreement with the non- 
interacting theory of Rcfs. 14,15 using a single fit param- 
eter, the dimensionless spin-orbit scattering rate A. The 
spin-orbit scattering time used in the fits was in order-of- 
magnitude agreement with previous measurements using 
weak localization.^^ The observation that the width of 
the distribution agreed well with theory after A has been 
chosen to fit the average was considered a success for 
random matrix theory. 

Significant deviations from the non-interacting theory 
can be expected for J/S > 0.2 only. Although this con- 
dition is not satisfied for the noble metals, the exchange 
interaction is strong enough to significantly affect the g 
factor distribution in most other metals, sec Fig. 9, where 
a list of values of J/S reported in the experimental and 
theoretical literature is given. The exchange interaction 
is particularly strong in Sc, V, Y, Nb, Rh, and, especially, 
in Pd (Pd is very close to the Stoner instability J/S = 1). 

For the collection of the data shown in Fig. 9, we used 
the fact that the ratio J/S is related to the Fermi liquid 
parameter Fq 



X = Xo-. 



7r 



pa ' 
^0 



(28) 



where m* is the effective electron mass, including band 
structure effects and interaction effects, m is the free 
electron mass, and xo is the Pauli susceptibility for free 

electrons,""^ 



Xo = 



(29) 



J/S = -F«, 



(27) 



The parameter 1/{1 + Fq) is also known as the Stoner 
enhancement parameter. 

The dimensionless spin-orbit scattering parameter A 
increases with element's nucleus charge Z. In experi- 
ments of Petta and Ralph^ strong spin-orbit scattering 
was found for Au nanoparticles of a few nm in diam- 
eter (A 0.1), whereas moderate spin-orbit scattering 
strengths (A ~ 1) were found in Cu and Ag nanoparti- 
cles of roughly equal size. From this, we conclude that 
moderate spin-orbit scattering strengths can be expected 
for nanoparticles with a value of Z around those for Cu 
or Ag. From Fig. 9 it can be seen that there are quite 
a few materials for which this is true and the criterion 
J/S > 0.2 is satisfied. 



V. CONCLUSION 

In this work we investigated the combined influence 
of electron-electron interactions and spin-orbit scattering 
on the g factors of metal nanoparticles. In the presence 
of electron-electron interactions, g factors must be at- 
tributed to (transitions between) many-electron states, 
instead of single-electron states. Many-electron states 
can have g factors larger than two, although these cannot 
be observed by tunneling spectroscopy because of "selec- 
tion rules" as long as spin-rotation symmetry is present. 
Spin-orbit scattering breaks spin-rotation symmetry and 
thus removes the "selection rules" . While this leads to a 
suppression of the g factors for large spin-orbit scattering 
rates, we find that g factors larger than two occur with 
significant probability if the spin-orbit scattering rate 
1/tso is moderate, TsqS < 1, where i5 is the mean spacing 
between single-electron energy levels in the nanoparticle. 
We have studied the g factor distribution quantitatively 
using random matrix theory and the universal interac- 
tion Hamiltonian.^'^ In addition to a confirmation of the 
scenario outlined above — occurrence of g factors larger 
than two — we found that interactions increase the width 
of the g factor distribution for t^oS 1 and that the g fac- 
tors probability distribution function is different for tran- 
sitions to the odd-electron ground state and odd-electron 
excited states. The enhanced fluctuations occurring for 
moderate spin-orbit scattering strengths (spin-orbit scat- 
tering rate 7so comparable to single-electron level spacing 



11 



1 H 








2 He 


3 Li 

0.23 


4 Be 

0.06- 
0.08 




5 B 


6 C 


7 N 


8 


9 F 


10 Ne 


llNa 

0.22 


12 Mg 

0.27- 
0.32 


13 Al 

-0.05- 
-0.2 


14 Si 


15 P 


16 S 


17 CI 


18 Ar 


19 K 

0.29- 
0.30 


20 Ca 

0.16 


21 Sc 

0.59- 
0.84 


22 Ti 

0.26- 
0.30 


23 V 

0.57- 
0.63 


24 Cr 


25 Mn 


26 Fe 


27 Co 


28 Ni 


29 Cu 

-0.03- 
0.09 


30 Zn 


31 Ga 


32 Ge 


33 As 


34 Se 


35 Br 


36 Kr 


37 Rb 

0.30 


38 Sr 

0.00- 
0.15 


39 Y 

0.53 


40 Zr 

0. 1 6- 
0.20 


41 Nb 

0.39- 
0.42 


42 Mo 

0. 1 3- 
0.19 


43 Tc 

0.24 


44 Rl 

0.25 


45 Rh 

0.41- 
0.44 


46 Pd 

0.84- 
0.89 


47 Ag 

0.02- 
0.07 


48 Cd 


49 In 


50 Sn 


51 Sb 


52 Te 


53 I 


54 Xe 


55 Cs 

0.18- 
0.36 


56 Ba 

0.22 


57 La 

0.52 


72 HI 
0.18- 
0.20 


73 Ta 

0.36 


74 W 
0.10 


75 Re 

0.15 


76 Os 

0.17 


77 Ir 

0.29 


78 Pt 

0.59- 
0.72 


79 Au 

0.06 


80 Hg 


81 Tl 


82 Pb 


83 Bi 


84 Po 


85 At 


86 Rn 


87 Fr 


88 Ra 


89 Ac 


104 


105 


106 


107 


108 


109 


110 


111 


112 


113 


114 


115 


116 


117 


118 



FIG. 9: Ratio J/5 of the exchange interaction constant and the mean level spacing. No value is listed for metals from Z = 24 

to 28 as they arc magnetically ordered as well as all lanthanides {Z = 58 . . . 71) except for Pm (Z = 61) for which no value 
was found in the literature. Also, for Z > 87 no data was found in the literature. For metals in the right-hand side of the 
periodic table (12th column and further) only data on their Pauli susceptibility in the liquid form arc available, sec Rcf. 26. 
For other metals, the data are taken from: Li, experiment;^^ Be, calculation of electronic structure;^^'^® Al, experiment^® 
(note the negative value of the exchange constant); K, experiments;^"'^^; Ca, Y, Tc, Ba, La, Ta, W, Re, Os, Ir, calculation 
of electronic structure;'''^ Sc, calculation of electronic structure;^^'^^ Ti, Zr, Hf, calculation of electronic structure;^'^'^'' V, fit 
of theory"^^ and experiment,'^'', and calculation of electronic structure;'^*'^^ Cu, Ag, experiment^"' and calculation of electronic 
structurc'^'^ (note the negative exchange constant in the experiment for Cu); Rb, Cs, experiment;'^" Sr, Nb, Mo, Rh, calculation 
of electronic structure;^'''^^ Pd, fit of theory^^ and experiment ;^'^ Pt, fit of theory"** and experiment,^® and calculation of 
electronic structure;^^ Au, calculation of electronic structure. 



i5) may bo typical of enhanced fluctuations of interaction 
matrix elements that are expected to occur at the break- 
down of random matrix theory (dimensionless conduc- 
tance ~ 1). 



There are ~ 20 metallic elements in the periodic ta- 
ble for which the electron-electron interactions are suffi- 
ciently strong that the phenomena described here can 
be measured. Existing measurements of g factors in 
nanoparticles have been made for Al and the noble met- 
als only; in these metals, interaction effects are weak. We 
hope that our findings stimulate experiments on other 
metals. 



In our calculations we have omitted the orbital contri- 
bution to the g factors. The orbital contribution arises 
from the fact that single-electron wavefunctions are com- 
plex if spin-orbit scattering is present, instead of real, so 
that they carry a finite current density. For the parame- 
ter regime we consider here, TsoS > 1, the single-electron 
wavefunctions are mostly real, and we expect the orbital 
contribution to the g factor to be smaller than the spin 
contribution. ^^'^^ We expect that taking into account a 
(small) orbital contribution causes a slight increase of the 
average g factor and a further broadening of the distri- 
bution. 
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APPENDIX A: g FACTORS IN THE WEAK 
SPIN-ORBIT SCATTERING LIMIT 



In this appendix, we calculate g factors to lowest order 
in the dimensionless spin-orbit scattering rate A. For a 
state that is twofold degenerate in the absence of spin- 
orbit scattering, spin-orbit scattering does not affect the 
spin contribution to g factors up to linear order in A.'*-'^ 
Our calculation addresses the case of a fourfold degener- 
ate level, and shows that its g factor is affected to zeroth 
order in A. This calculation shows explicitly that the 
limit A ^ is singular. 

Let us start by writing the spin-orbit Hamiltonian in 
terms of the basis of single-electron eigenstates of the 
Hamiltonian Hqoe without spin-orbit scattering, see Eq. 
(4) above, 



iX 
2s/N 



+ (Al + iA2)f,.i^l^i^,^ 
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FIG. 10: An S = 3/2 state and the three relevant types of 
excited states; (a)the four-fold degenerate 5 = 3/2 ground 
state; (b) two-fold degenerate 5 = 1/2 excited state; (c)four- 
fold degenerate S = 3/2 excited state; (d) four-fold degener- 
ate 5=1/2 excited state, note that this state is entangled 
see Eq. (A6) 



(Al) 



Here the Greek indices fj, and u refer to the eigenvalues 
of Hgok, and not to the eigenvalues of the total single- 
electron Hamiltonian Hqqe + Hso as in Sec. 11. 

In order to study the effect of spin-orbit scattering on 
the g factor of a many-electron eigenstate of Hqoe with 
spin S = 3/2, one needs to calculate matrix elements 
between the four members of the quadruplet. Labeling 
the four members of the quadruplet by the z component 
of the spin, Sz = p — 5/2, p = 1,2,3,4, these matrix 
elements can be arranged in a 4 x 4 matrix V of the form 



V 




b c 

-a + d c 
-a + d -b 

c* -b* -a-d, 



(A2) 



with a and d real numbers and b and c complex numbers. 
The specific form of (A2) follows from time-reversal sym- 
metry and guarantees that that the eigenvalues of V are 
double degenerate, in accordance with Kramers' theorem. 

One quickly verifies that all matrix elements of V are 
zero to first order in Hso ■ This is the consequence of the 
fact that the matrices Aj, j = 1, 2, 3, in Eq. (Al) are an- 
tisymmetric, so that the spin-orbit interaction does not 
mix the states with opposite spin belonging to the same 
energy level. In this situation one has to calculate ele- 
ments of V to second order in i?so-^^ Denoting the many- 
electron states with roman indices, the matrix elements 
between the many-electron states n and n' (both taken 
from the same quadruplet) are given by 



E 



(-^so ) Tim (-^so ) 7T 



(A3) 



with m is summed over all many-electron states with 



and E^^^ 



Em 7^ 

electron energies. 



and the corresponding many- 



The quadruplet state is represented schematically in 
Fig. 10a. For the calculation of the splitting of a quadru- 
plet, it is enough to consider states m of the form indi- 
cated in Figs. Fig. 10b, c, and d. These are: a twofold 
degenerate S = 1/2 state (Fig. 10b), a fourfold degener- 
ate 5 = 3/2 state (Fig. 10c) and a fourfold degenerate 
S = 1/2 state (Fig. lOd). There exist two variants of the 
states shown in Fig. 10c and d, depending on whether an 
empty single-electron level is filled above the Fermi level, 
or a hole is created below the Fermi level. The former 
case is shown in the figure. Since the spin-orbit inter- 
action is a one-particle operator, only the states of the 
form shown in Fig. 10 which differ by not more than one 
electron-hole excitation are important. 

Virtual excitations to two-fold degenerate 5 = 1/2 
state. In this case, the transition from the state n to 
m involves a transition of an electron from the singly- 
occupied level jjL to the already (singly) occupied level u, 
see Fig. 10a and b. Representing the members of the spin 
S = 3/2 quadruplet as \'i/2,S^) with = -3/2, -1/2, 
1/2, 3/2, and the members of the spin 5=1/2 doublet 
as 11/2,52) with Sz = 1/2, —1/2, we find the following 
matrix elements of the spin-orbit Hamiltonian Hso, 

(3/2,+3/2|if,o|l/2,+l/2) = -^(Ai-^a)^,, 

(3/2,+3/2|if,o|l/2,-l/2) = 0, 

(3/2,+l/2|iJ,o|l/2,+l/2) = 

(3/2,+l/2|if,o|l/2,-l/2) = 

(3/2,-l/2|i?,o|l/2,+l/2) = 

(3/2,-l/2|if,o|l/2,-l/2) = 
(3/2,-3/2|if,o|l/2,+l/2) = 
(3/2,-3/2|/fso|l/2,-l/2) = 

Further, in this case the energy difference 

- =e^-e^ + 3 J. (A5) 

Substituting these matrix elements and the energy dif- 
ference into Eq. (A3) we find the following contributions 
to the elements of the matrix V of Eq. (A2): 




4N{e, -e^ + 3J)\/3 



4A^(£, . 



3J)y3 



X\A,)l^+X\A,)l^-2\\A3)l^ 



127V (e^ -e^ + 3 J) 

(We have not listed the value of a in the matrix V of Eq. 

(A2) since this coefficient does not contribute to the g 
factor and splitting of the 5 = 3/2 quadruplet.) 
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Virf,ual excitations to four-fold degenerate 5 = 3/2 
state. These excitations involve a transition of an elec- 
tron from a singly occupied level fj, to an unoccupied level 
u, see Fig. 10, or the transition of an electron from a 
doubly occupied level /x to a singly occupied level p. Cal- 
culating the various matrix elements as before, we find 
after somewhat cumbersome algebra 



C2,au = — 



2,1X1/ 



\^{Ai-iA2)^,„{A?,)fj,u 
6N{e, - e^)V3 

6N{s^ - £^)^/3 

XHA^)l^ + X''{A2)l^-2X\As)l^ 



Virtual excitations to four-fold degenerate S — 1/2 
state. As in the previous case, these excitations involve a 
transition of an electron from a singly occupied level n to 
an unoccupied level see Fig. 10, or the transition of an 
electron from a doubly occupied level /x to a singly occu- 
pied level u. The four S = 1/2 are labeled by = ±1/2 
and by an additional degeneracy parameter q = ±1, 



|l/2,+l/2,g) 



TTi) + 1 Tine'^S^ 



V3 ( 

+ I iTT>e-^''* 

1 ( 



11/2,-1/2,(7) = ^(|UT) + |m)e'^ 

+ I Tii)e-''^) . (A6) 
Performing the calculations as before, we find 
X'^{Ai-iA2)^,„{Aa)f,^ 



b3,u.v = 



d3,ij,v 



3J)V3 



X\A^-iA2)% 
6N{s^ -e^ + 3J)\/3 
XHA^)l, + X^A2)l,-2XHAs)l^ 



3 J) 



Denoting the set of doubly occupied single-electron lev- 
els by "0" , the set of singly-occupied single-electron levels 
by "1" and the set of unoccupied single-electron levels by 
"2" , we then sum over all virtual excitations and find 

(A7) 

and similar expressions for the coefficients c and d in Eq. 
(A2). 

The matrix (A2) can be diagonalized for all values of 

the parameters a, 6, c, and d, and the corresponding g- 
factors can be found exactly. After diagonalization we 



find that the quadruplet is split into two doublets with 
energy separation 



{AEf =4d2+4|6|2 + 4|c|' 



(A8) 



In the 4x4 matrix notation of Eq. (A2), the Zeeman 
Hamiltonian reads 



Hz 



\ 



HbH 



(A9) 



3hbH J 



Lifting the degeneracy of the two doublets by the Zeeman 
energy, we find g factors 



5 = 24 



l3\b\^ + {2d±AEy 
d2 |6|2 + |c|2 



(AlO) 



where the ± sign refers to the two doublets. This result 
confirms the assertion made earlier, that spin-orbit scat- 
tering affects the g factors of the S — 3/2 states to zeroth 
order in the spin-orbit scattering rate A. Of course, for 
small A the energy splitting AE is small as well, and the 
g factor of Eq. (AlO) can be observed for magnetic fields 
such that HbH ^ AE only, which limits the practical 
observability of the g factor (AlO) for very small spin- 
orbit scattering rates A <C 1. 

In the remainder of this appendix we investigate Eq. 
(AlO) for the special case that one state is very close to 
the 5 = 3/2 state of interest and virtual excitations to 
that state dominate the spin-orbit matrix V in Eq. (A2). 
One important example is the case when the ground state 
has spin S = 3/2, which is expected with small proba- 
bility for J/S > 0.3, sec Fig. 11. Indeed, the energy dif- 
ference between the 5 = 3/2 state and the lowest- lying 
5=1/2 state is £\+2 — £\ — 3 J, X being the index of 
the lowest singly occupied level in the 5 = 3/2 state. 
Typically, this energy difference is small, given the small 
likelihood of it being positive. 

Substituting the general expressions for the coefficients 
b, c, and d into Eq. (AlO) and noting that one energy 
denominator is much smaller than all others, we find that 
an 5 = 3/2 quadruplet splits into two doublets with g 
factors 



9 = 



9 = 




Aj + Al + Ai 



(All) 



where we have omitted the indices referring to the levels 

fi and v involved as we deal with one excited state only. 
If the quadruplet state is the ground state, g corresponds 
to the lower lying doublet. 

In the special case of an 5 = 3/2 groimd state, Eq. 
(All) shows that the (^-factor takes values in the interval 
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"^S?" has spin 1/2, wc expect the true average ground state g 

factor to be approximately equal to 



<g > 
<g> 



(a) 



1-P. 



3/2 



(b) 



FIG. 11: Two different possibilities: (a) S = 3/2 state is the 
ground state with the probability (b) S = 1/2 state is 

the ground state with the probability 1 — ^3/2- Note that due 
to the sign change in the energy denominator, see Eq. A3 the 
average g-factors < g > and < g' > are inverted. 



[2, 4] only. One has 5 = 2 only if Ai ^ A2 = 0. The 
corresponding eigenstates are |3/2,±l/2). In the oppo- 
site case A3 = 0, one has = 4 in the ground state and 
the corresponding eigenstates are (\/3/2)|3/2, ±3/2) — 
(e2*'^/2)|3/2,Tl/2), where e"^"^ = {A2+iAi)/\A2+iAi\. 
In fact, one can find the entire probability distribution 
of g in this case, using the fact that Ai, Ai, and A3 are 
taken from identical and independent Gaussian distribu- 
tions, 



Pis) 



1 



2\/3^16-52" 



(A12) 



The average g'- factor for an S = 3/2 ground state with a 
nearby S =1/2 state is then 



{9) 



4 fj: 
3.418. 



arcsm - 
2 




(A13) 



For J = 0.6(5, the probability to find a ground state with 
5 = 3/2 is P3/2 « 0.38 (see Ref. 21 or Fig. 8). Since the g 
factor is unaffected to first order in A if the ground state 



{go) « 3.4I8P3/2 + 2(1 - P3/2) « 2.54. 



(A14) 



This value is very close to that found in the numerical 
simulations of the ground state 5-factor w 2.5, see Fig. 4. 

The distribution of g' can be found from Eqs. (All) 
and (A12). One finds the particularly simple result 



< 5' < 6. 



(A15) 



The 

r,' 



average {g') = 3. In the hmiting cases 
g' = 6 and g' = the doubly-degenerate eigen- 
states have the form |3/2,±3/2) and (l/2)|3/2, ±3/2) - 
(y3/2)e2*'^|3/2,=Fl/2), respectively. We can use these 
results to calculate the average ^-factors of the lowest 
excited states. For the first excited state 

{91) = {9')P3/2 + {9') (1 - P3/2) = {9') = 3, (A16) 
where we assumed that without spin-orbit scattering the 
ground state has spin 5 = 3/2 or the the first excited 
state has spin S = 3/2 and is slightly below a.n S = 1/2 
doublet, see Fig. 11. (Hence, we neglect the possibilities 
that the first excited state has spin 5 = 1/2 or that 
the first excited state has spin S = 3/2 and is far away 
from the next S' = 1/2 state. Our approximation should 
slightly overestimate {gi)-) The simulations give {gi) w 
3.00, see Fig. 4. With the same approximations, we find 
that the average g factor of the second excited many- 
electron state is 



(52) = 2P3/2 + {9) (1 - ^^3/2) « 2.87. 



(A17) 



The result of the simulation is 2.78, in good agreement 

with the estimate (A17). Note that the above nontriv- 
ial distributions of (^-factors can be observed in a small 
magnetic field only {hbH <C For larger fields the 

g-factors are the same as in the absence of spin-orbit cou- 
pling. 
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